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ABSTRACT 

Using a parameterised function for the mass loss at the base of the post-shock region, 
we have constructed a formulation for magnetically confined accretion flows which 
avoids singularities, such as the infinity in density, at the base associated with all 
previous formulations. With the further inclusion of a term allowing for the heat 
input into the base from the accreting white dwarf we are able also to obtain the 
hydrodynamic variables to match the conditions in the stellar atmosphere. (We do 
not, however, carry out a mutually consistent analysis for the match). Changes to the 
emitted X-ray spectra are negligible unless the thickness of mass leakage region at the 
base approaches or exceeds one percent of the height of the post-shock region. In this 
case the predicted spectra from higher- mass white dwarfs will be harder, and fits to 
X-ray data will predict lower white-dwarf masses than previous formulations. 
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1 INTRODUCTION 

. In magnetic cataclysmic variables (mCVs) the accreting material from the secondary star is entrained onto magnetic field 
lines and accretes onto the surface of the white dwarf. There it forms a standoff shock, followed by a hot post-shock region 
of plasma settling onto the white dwarf surface as it cools, principally by bremsstrahlung X-ray radiation, and by optical/IR 
cyclotron radiation if the magnetic field is sufficiently strong (see Cropper 1990 and Warner 1995 for reviews of mCVs). 

There have been a number of studies of the post-shock accretion flow (e.g. Aizu 1973; Langer, Chanmugam & Shaviv 
1981; Chevalier & Imamura 1982; Kylafis & Lamb 1982; Chanmugam, Langer & Shaviv 1985; Imamura etal. 1987; Wu 1994; 
Wu, Chanmugam & Shaviv 1994; Woelk & Beuermann 1996; Kocabiyik 1997; Cropper etal. 1999; Saxton & Wu 1999), 
investigating different aspects of the region, both analytically and numerically (see Wu 2000 for a review). Despite the fact 
that substantial progress has been made, the existing studies generally assume that the velocity and the temperature drop to 
zero at the base of the accretion column. The assumption of zero flow velocity and the strict requirement of mass continuity 
along the field lines immediately implies that the matter density at the base of the accretion column must reach infinity. 
Recent work by Cropper, Wu & Ramsay (2000) has highlighted the consequences of an infinite density at the base of the 
accretion column: the X-rays from the shock-heated region is emitted mostly from its base where the density is rising steeply. 
This is true at energies up to approximately that of the shock temperature, which is in the ~ 10 — 60 keV range — higher 
than most imaging or CCD-based X-ray instruments. The preponderance of emission from the base of the accretion column 
has made a more in-depth understanding of the physical conditions in shock-heated emission regions essential, not only from 
the theoretical point of view but also for modelling and extracting information from the observed X-ray spectra. 

Several aspects have to be considered in order to derive a self-consistent formulation of the accretion flow, especially at 
the boundary layer where the hydrodynamic flow merges into a hydrostatic white-dwarf atmosphere. In the stationary-state 
case, the hydrodynamic variables at the base of the accretion column should match smoothly to the corresponding variables 
at the white-dwarf atmosphere. In addition, the energy deposition due to the heat flux emerging from below the white-dwarf 
atmosphere should be considered. Moreover, at some point in the post-shock region, material may become sufficiently cold 
that a fraction becomes neutral and cannot be efficiently confined by the magnetic field. In the time-dependent situation, it 
is necessary to consider the response of the white-dwarf atmosphere to changes in the local mass-accretion rate and its effects 
on the stability of the accretion flow. 

In this paper we investigate the stationary-state accretion onto magnetic white dwarfs, with emphasis on the post-shock 
flow at the boundary layer above the white-dwarf atmosphere. We consider a modification to the conventional hydrodynamic 
formulation so that matter diffusion across the field lines is allowed in a thin region at the bottom of the accretion column. 



2 Wu & Cropper 



In addition, the heat flux from the white-dwarf atmosphere into the accretion column is considered. As it is an exploratory 
study, we neglect the complications such as geometrical, gravity and two-fluid effects. We demonstrate that a closed-form 
solution as that described in Aizu (1973) and Wu etal. (1994) (see also Chevalier & Imamura 1982) can be obtained for this 
formulation. Moreover, the infinities and discontinuities at the base, which lead to unphysical observational consequences, are 
eliminated. We calculate the emission from the modified post-shock region and compare their spectral properties with those 
of the conventional post-shock region. 



2 STATIONARY-STATE ACCRETION FLOW 
2.1 Modified hydrodynamic formulation 

We assume that the flow of the ionised accreting material follows the magnetic field lines and that the field lines are per- 
pendicular to the white-dwarf surface, so the flow is (quasi-)one-dimensional. We neglect gravity effects, which are important 
when the shock height x s is not negligible in comparison with the white-dwarf radius R w (see Cropper etal. 1999). We also 
assume that only the ionised matter is strictly confined by the magnetic field, while the cold neutral atomic matter can "leak" 
out of the accretion column. The accretion matter obeys the ideal gas law, i.e. P = 2pkT/m p and 7 = 5/3, where P is the 
gas pressure, T the gas temperature, k the Boltzmann constant, m p the proton mass and 7 the adiabatic index of the gas. 
At the lower boundary of the accretion column the temperature gradient and the bulk-flow velocity are zero, and the cooling 
rate equals the heating rate. 

In our formulation the stationary state mass continuity, momentum and energy equations are 
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where v is the flow velocity, and p the density. E is the sink term specifying the rate of mass loss out from the accretion 
column; A c is the effective cooling function, and Ah is the heating function. In the construction of the formulation, we have 
implicitly assumed that the neutral matter, which is not confined by the magnetic field, transfers its momentum and energy 
to the ionised matter before it leaves the accretion column. Under this assumption momentum is conserved along each flow 
line of the ionised material and energy is dissipated only by emitting radiation. In reality the neutral matter may carry away 
substantial amount of momentum and energy, and hence the strict condition for conservation of momentum and energy along 
the field lines needs to be relaxed. This case will be discussed elsewhere (Wu & Cropper, in preparation) . 

The accretion shock is assumed to be strong and adiabatic, with the immediate post-shock velocity equal to a quarter of 
the pre-shock velocity, which is taken to be the free-fall velocity at the white-dwarf surface, i.e. v — —vg/A at x = x a . (Here 
and hereafter, the subscript "s" denotes variables at the shock.) The post-shock pressure is three-quarter of the pre-shock 
ram pressure of the inflowing gas, i.e. P B = 3p a v|/4, where p a , the density of the pre-shock accretion matter, is related to the 
specific accretion rate (per unit area) m by p a = m/vg. 

We define the dimensionless variables £ = x/x B , £ = p/p a , t = —v/vs and vo = P/p^Vg, and substitute them into the 
mass, momentum and energy equations. In terms of these variables, the hydrodynamic equations are 
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where E = (x B / p a vg)T,, and A = (x s / p a fg)(A c — A h ). At the white-dwarf surface (£ = 0), the velocity r = 0; at the shock 
(£ — 1), t — 1/4 and the pressure vj — 3/4. 

The sink term above, E, is not a direct observable, but its integration over the distance along the flow is the mass loss 
rate per unit area. It is therefore more appropriate to consider the integration of the sink term over the distance and express 
it in terms of the flow velocity, the independent variable that we use in solving the hydrodynamic equations. 

In terms of the integration of the sink term over the distance we can define a dimensionless variable 
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At the white-dwarf surface <r(0) = 0, and at the shock <r(l/4) = 1. The variable <t(t) clearly satisfies the mass continuity 
equation and it is in fact the dimensionless specific mass accretion rate in the presence of matter leakage. 
The dimensionless density can now be obtained by solving the mass continuity equation, and it is 
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The dimensionless pressure, obtained by integrating the momentum equation, is 

ot(t) = 1 — T <t(t) . 

Combining equations (8) and (9) with the energy equation yields 
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The above equation can be integrated to obtain a closed-form expression for £ (t) 
i(r) = 1 - J'" dr'l 7 - (7 + 1) r' a - \ (7 + 1) r' 2 A a 

We follow Wu (1994) and consider the cooling term A c (in equation [3]) in a composite-functional form 
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The parameter e s specifies the ratio of the bremsstrahlung cooling time-scale to the cyclotron cooling time-scale at the shock. 
The value of the constant A is 3.9 x 10 16 in c.g.s. units for pure hydrogen plasmas (see Rybicki & Lightman 1979), and a ^ 2 
and P ^ 3.85 for parameters appropriate for mCVs. The dimensionless cooling/heating term A is then 
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With the cooling function is defined, the flow velocity, and hence the other hydrodynamic variables, can be derived after 
equation (11) is inverted. 

If we set 7 = 5/3, o(t) — 1 and Ah = 0, equation (11) becomes 
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which is identical to equation (3) in Wu (1994). This limiting case corresponds to "no leaking" from the accretion column 
and a hard, stationary white-dwarf surface with no emergent energy flux from below. If we further set e s = 0, the Aizu (1973) 
solution is recovered. 



2.2 The sink term 

When the accreting material is cooled sufficiently, recombination occurs and neutral atoms are formed. While the charged 
ionised matter is magnetically confined, neutral matter can diffuse across the field lines. As the temperature is lower near the 
base of the accretion column, we expect matter leakage to be efficient. 

In the hydrodynamic formulation that we consider, the temperature and the flow velocity are not independent. The 
temperature is a monotonic function of the velocity, and in the limit of no leakage it is linearly proportional to the velocity 
at the base of the accretion column. The sink (which depends on the hydrodynamic variables) can therefore be chosen to be 
an explicit function of the flow velocity only. We assume that it is characterised by a critical velocity r m , below which matter 
leakage becomes efficient. This is equivalent to assuming a critical temperature at which recombination becomes sufficiently 
efficient to allow significant amount of neutral matter to form and to diffuse out of the accretion column. 

It is more convenient to consider <t(t) (the specific mass- accretion rate per unit area remaining in the post-shock flow at 
height £ after leakage) than the leakage at any height E(£) itself. Suppose <t(t) can be expanded into a series of orthogonal 
bases e~ nT//Tm , such that <t(t) = ^2 a, n e^" T ^ Tia , where a n is the coefficient of the n-th power term. If only the first two leading 
terms are important, as an approximation we may consider 

a(r) = p (l-ve-^), (15) 

where p and v are constants to be determined. In this study we treat r m as a parameter. When appropriate atomic and 
magneto-hydrodynamic processes are considered, an explicit expression of it in terms of the other system parameters can be 
obtained. 

The functional form of <t(t) above must satisfy the hydrostatic-equilibrium condition at the white-dwarf surface (£ = 0). 
This requires the mass flux o(j) equal to zero at the white-dwarf surface. Thus, the constant v = 1. At the shock (£ = 1), the 
mass flux o(j) = 1 by definition (equation 7), implying a normalisation constant p = (1 — e -1 / 4rm ) -1 . In the limit of r m <JC 1, 
the velocity-derivative of <t(t) is 
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Figure 1. Specific mass flux <x(£) (mass transfer rate per unit area) as a function of height £ within the post-shock region for T m 
0.001, 0.01 and 0.05 (curves a, b and c respectively). 



at the white-dwarf surface, and it is 
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at the shock. Matter leakage therefore occurs only at the very bottom of the accretion column if r m is sufficiently small, and 
mass is practically conserved along the flow in most of the post-shock region. For r m = 0.01, the specific mass flux falls to 
less than 1/2 of its initial value only when the normalised height £ £ 0.001; and for r m = 0.001, when ^3x 10~ 6 (Fig. 1). 



2.3 Energy flux from the white-dwarf atmosphere 

We assume that the heating due to the energy flux from below the white-dwarf surface is important only in a thin region (with 
a thickness At < r m ) at the bottom of the accretion column. Suppose ttSv,(0) is the radiation flux emerging normally from 
the white-dwarf surface. The radiation is absorbed by the accreting matter and is attenuated. At the height £, the radiation 
flux is given by 



ttS w (£) = ttS w (0) exp 
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where is the absorption coefficient. The rate of energy deposited per unit volume is therefore 
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The absorption coefficient k depends only on the local density and temperature, which are approximately constant when 
t < Tm (see next section). In local thermal equilibrium, the absorption coefficient is determined by the emissivity, which is 
the effective cooling function in our formulation, and the Planck function B (Kirchoff's Law). As the emergent flux from the 
white dwarf is a black-body flux, i.e. S , w (0) = 13(0), we have 
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We have defined the lower boundary of the accretion column as the location where the heating rate equals the cooling 
rate and the flow velocity is zero. This requires Ah(r) = A c (r) and r = at £ = 0. Suppose the attenuation factor can be 
expressed in terms of the dimensionless velocity r, such that 

KM 

(22) 



where A is a constant to be determined. Then, the heating function in the boundary layer is 
A h (r) « A c (0) exp 



A 

By expanding r into Taylor series at the boundary £ = 0, we obtain 
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Smooth merging of the flow into a hydrostatic white-dwarf atmosphere requires the first derivative of the velocity to be zero. 
By combining equations (10), (23) and (24) and keeping only the first non- vanishing term, we obtain 
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We can now define a characteristic critical velocity r h , below which heating by the flux from below the white-dwarf 
atmosphere is important and beyond which heating effects are negligible. Th is given by 
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The heating function 
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is asymptotically zero for r ^> Th. As heating is unimportant for r > r m > rh, we can simply use the heating function above 
throughout the whole post-shock region, despite the fact it is derived from considering heating of the accreting material within 
the boundary layer. 



3 PROPERTIES OF THE HYDRODYNAMIC VARIABLES AT THE LOWER BOUNDARY 
3.1 Hydrodynamic variables 

The matter density at the base of the "leaky" accretion column that we consider is 
C(0) = lim H (l - e" T/Tm ) 

T m 



(28) 



which is clearly finite. The pressure is 
B7(0) = lim h - rn (l - e- T/ Ml 

= 1, (29) 
the same as that in the conventional formulation. The dimensionless temperature is defined as 9 = T/T B , and hence we have 



At the white-dwarf surface, the temperature is 
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which is also finite. 

In the limit of r m — » 0, we have er(f) — » 1 and da/dr — » 0. Moreover, ("(0) — » co, ^(0) — + 0, and tu(0) = 1, i.e. the cold 
"stationary-wall" boundary condition is recovered. 
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3.2 Gradients of the hydrodynamic variables 



Not only are the hydrodynamic variables of the conventional and our formulations are different at the lower boundary, but 
they also have very different gradients. In the conventional formulation, the density gradient is given by 
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As dr/d£ oc r 3//2 for small r (from equation (2) in Wu 1994), the density gradient is proportional to r 7 ^ 2 . At £ = 0, r = 0; 
hence the density gradient is infinite. The pressure gradient, which is 
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and the temperature gradient, which is 
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are also infinite at £ = 0; 

For a "leaky" accretion column allowing an energy flux emerging from below, the velocity gradient at the base is 
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The velocity gradient is zero at £ = 0, provided that A c (0) = Ah(0). This is in fact the condition of radiative equilibrium in a 
hydrostatic white-dwarf atmosphere. 

Differentiating equation (8) with respect to £ yields the density gradient: 
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By substituting the expression of <t(t) given in equation (15) into the equation above, we obtain 
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The power series C n converges, with a value of the order of unity when t — ► (Appendix A.l). The density gradient is 
therefore finite and equals zero at £ = 0. 
The pressure gradient is 
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It is also finite and equals zero at £ = 0. The temperature gradient is 
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As dzu/dr and d(/dr are both finite and the density £ is non-zero at £ = 0, the temperature gradient is zero. 



4 ACCRETION LUMINOSITY 

In the conventional formulation, energy conservation holds along the field (flow) line, so that the accretion energy must be 
radiated away before the accreting matter can settle down onto the white-dwarf surface. The accretion luminosity (normalised 
to p a w|) along the flow line down to the height £ is 
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Figure 2. Density f (£) as a function of height £ for T m = 0.001, 0.01 and 0.05 (curves a, b and c respectively). The density dependence 
of the conventional model (Aizu 1973) is also shown for comparison (dotted line). 
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For 7 = 5/3, L acc (0, 1) = 1/2, consistent with the fact that the total accretion energy liberated is independent of the nature 
of the cooling process. 

For the "leaky" accretion column considered here, energy conservation also holds strictly along the flow lines. The 
bolometric accretion luminosity along a field line down to £ is 
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(Appendix A. 2), where a is as defined in equation (15), and the total bolometric accretion luminosity is 
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For 7 = 5/3 it is equal to 1/2, consistent with the requirement of strict energy conservation along a field line. 
The bolometric accretion luminosity of the boundary layer, where r < r m , is 
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where £ m is the height at which r = r m . If the boundary layer is thin (i.e. r m 1), the bolometric luminosity ~ 7T m , which 
is insignificant in comparison to that of the rest of the post-shock region. 



5 DISCUSSION 
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Figure 3. Density (left) and temperature (right) as a function of height £ for heating scale height = 0.001 and 0.0001 (curves a and 
b respectively). The leakage parameter T m is fixed to be 0.01. The conventional model is shown (in dotted line) for comparision. 



5.1 Structure of the post-shock region 

Figure 2 shows the density profile at the base of the post-shock region for "leaky" accretion columns with different values 
of Tm by comparison with the conventional formulation (Aizu 1973), in which the lower boundary is a cold, stationary wall. 
As noted above, the density £ and its gradient in the conventional formulation approach infinity at the base of the accretion 
column, whereas the density in our formulation has a finite value, ^i/r m . For r m = 0.05, 0.01 and 0.001, £(0) « 20, 100 and 
1000 respectively. In all cases of r m , the density at the "leaky" base falls significantly below that of the base in the conventional 
formulation for £ less than £ m , the height at which r = r m . 

In calculating the density profiles in Figure 2, we have not included the effect due to heating from the white dwarf, i.e. 
we have considered Th = 0. We now show this effect in Figure 3, where we have fixed the critical velocity r m to 0.01, and 
calculated the dimensionless density and temperature for two different boundary layers, with r h = 0.001 and 0.0001. They 
both fulfill the constraint that r h < r m . As shown, the terminated temperature depends mainly on r m . The heat input from 
the white dwarf modifies only the asymptotic properties of the base temperature — a larger Th will result in a slightly thicker 
"isothermal" layer at the base, thus reducing the temperature upstream and causing the density to reach the terminated 
plateau value further upstream. 

Figure 4 shows the dimensionless local bremsstrahlung emissivity jbi (= C^O 1 ^ 2 ) as a function of the height £ in the post- 
shock region. The emissivity reaches a terminated value determined by r m , whereas that from the conventional formulation 
tends to infinity. Again the bremsstrahlung emissivity is modified slightly by the input of heat from the white dwarf when 
Th > with increased emission for larger values of Th. This is in fact the consequence of an implicit assumption that the 
white-dwarf energy flux is much weaker than energy flux released by accretion material. 



5.2 Optical depths 

Using the prescription for density and temperature above we can easily show that the optical depth due to free-free absorption 
or electron scattering in the vertical upstream direction is negligible when the photon energies are above 0.05 keV, for typical 
mass-transfer rates and white-dwarf masses. However, by assigning the parameter Th for the heating at the base of the region 
we implicitly assume a large opacity to maintain a black-body spectrum for the radiative flux. As the temperature in the 
region of our interest should be sufficiently low such that recombination can occur to form neutral atoms, there are additional 
sources of opacity, such as the bound-free and bound-bound opacities, to maintain a large optical depth in this geometrically 
thin layer at the base. 

The bound-free and bound-bound emission processes are, however, not considered explicitly in constructing the heating 
and cooling functions in our formulation. Instead, we consider only a parametric heating function Ah, which is in terms of 
the accretion parameters and the energy flux from below the white-dwarf surface. The formulation that we present in this 
work is therefore not a fully self-consistent formulation with explicitly consideration of the microscopic atomic physics (see 
van Teeseling, Heise & Paerels 1994). 

Also, we have not considered an explicit treatment of the ionisation equilibrium. There might be situations that the gas at 
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Figure 4. Bolometric bremsstrahlung emissivity jb T (= f 2 ^ 1 / 2 ) as a function of height £. Curve a corresponds to the case of r m = 
0.01 and Th = 0.001; curve b to the case of r m = 0.01 and = 0. The bremsstrahlung emissivity of the conventional model is shown (in 
dotted line) for comparison. 



the base is photo-ionised by the X-ray from the shock above. The increase in the degree of ionisation will cause an decrease in 
the efficiency of matter diffusion across the field lines, and hence out of the accretion column to the white-dwarf atmosphere. 
For typical mCV parameters, say white-dwarf mass M wd = 0.7 M Q , specific mass-accretion rate m ~1 g cm _2 s _1 , X-ray 
luminosity L x of ~ 10 32 erg s _1 , shock temperature T B ~10 keV and shock height i s ~3x 10 7 cm, the ionisation parameter 
S (= L x /nr 2 , where n is the number density and r is the distance to the X-ray source) has values ,$,0.01. (Here, we have 
assumed r m ~ 5.0 x 10~ 5 in the calculation of H.) The value of H is significantly less than 20, the critical value at which 
hydrogen becomes completely ioinised (see Models 3 and 4 in Kallman & McCray 1982). Hence, the formulation presented 
above is in general applicable to the accretion column of mCVs. Such a parametric prescription allows to derive a analytic 
formulation to describe the accretion column and the boundary layer in a simple and unified manner. 

It is worth noting that the surface at which the total scattering and absorption optical depth is unity is in effect the 
observable lower boundary of the post-shock region. This surface does not necessarily coincide with the base as that defined 
by the flow hydrodynamics, the zero-velocity surface. The emission from this unit-optical-depth surface will be observed as 
soft X-rays with a blackbody-like spectrum characterised by a temperature similar to that of the isothermal layer (due to 

Th < T m ). 



5.3 X-ray spectral properties 

With the density and temperature specified within the post-shock region we can calculate the emission at different heights 
and construct integrated spectra. We consider the case in which bremsstrahlung radiation is the only cooling process and 
ignore cyclotron cooling, so to illustrate the differences between between it and the conventional post-shock X-ray emission 
regions resulting from the Aizu (1973) formulation. 

Figure 5 illustrates the differences in the emitted spectrum with respect to the conventional formulation for I.O-M0 and 
0.5-Mq white dwarfs. The variable tpbf(E) is the ratio of the bremsstrahlung emissivity for the "leaky" post-shock region 
jbi(E; r m ) to that for the conventional region jbi(E; 0), where E is the photon energy. When r m is relatively large (~ 0.01), a 
significant fraction of accreting material is lost from the post-shock region at large height (see Fig. 1) . Because we require all 
of the accretion energy to be liberated before the matter diffuses out of the post-shock region, the accretion energy is radiated 
at higher energies and the 1-20 keV region of the resulting spectrum is enhanced with respect to that of the spectrum from 
the conventional post-shock region. 

For 1.0-Mq white dwarfs the spectrum below ~ 5 keV is harder than that from the conventional post-shock region, while 
above this value it is softer. This effect is enhanced if significant heat flux is permitted from the white dwarf (larger values 
of Th). For lower values of r m , the temperatures are already low at the base of the region, and the main effect is that of 
the increased heating from the white dwarf at energies below 1 keV. The spectrum in this case is softer than that from the 
conventional formulation, but only at low energies. When both r m and r h tend to zero, the spectrumretains the conventional 
form. 
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Figure 5. The ratio of total brcmstrahlung cmissivity of the "leaky" accretion column to that of the conventional accretion column ip^, T 
as a function of photon energy E for a 1.0-Mq white dwarf (left) and a 0.5-Mq white dwarf (right). The leakage and heating parameters 
are: (r m ,T h ) = (0.01,0), (0.01,0.001), (0.001,0.0001) and (0.0001,0.00001) for cuves a, b, c and d respectively. The specific mass flux is 1 
g s -1 cm -2 at the shock (£ = 1). 

5.4 Observational consequences 

From Figure 5 it is clear that it is possible to obtain either a harder or softer spectrum than that from the conventional 
post-shock region, depending on leakage parameter r m and heating parameter Th, and on the energy range considered. When 
fitting to X-ray spectral data, a harder model spectrum will result in a lower-mass determination for the white dwarf. If r m is 
significant, then for typical CCD or proportional counter detectors operating in the 0.2-20 keV range, fitting spectra from the 
"leaky" formulation will generally therefore result in lower masses for more massive white dwarfs than those obtained using 
the Aizu (1973) or Wu ct al. (1994) formulations. In the case of lower-mass white dwarfs the effect is to mimic the change in 
spectral slope resulting from increased cyclotron cooling from a strong magnetic field. 

It should be noted that because they are parameterisations, the treatments in §2.2 and §2.3 do not provide a prescription 
for definite values of T m and Th based on atomic physics and magneto-hydrodynamics. However, if we simply assume certain 
values for the effective temperature T e g of the white-dwarf atmosphere, then we can obtain an estimate for r h by equation 
(26) after the specific accretion rates (or the specific pre-shock mass flux) m and the white-dwarf mass M w are specified. If we 
further assume that the temperature T(0) at the boundary (i.e. where r and its derivative is zero) as the unperturbed effective 
temperature of the white-dwarf atmosphere, then we may also obtain an estimate for r m using equation (31). In reality, we 
expect T(0) to be larger because of radiative heating by the X-rays from the post-shock region above and T B to be less than 
that determined by the free-fall velocity at the shock surface (see Cropper etal. 1999). Therefore, we have underestimated r m 
here, and the value that we obtained can be considered as a lower limit to r m - Nevertheless, under this simple approximation, 
we have found that the assumption Th < r m is acceptable, provided that T e g is not very significantly larger than 10 4 K. 
Moreover, the assumption that r m <C 1 is in general valid. This is shown in Figure 6. 

Which values of r m can be reached is therefore unclear pending further investigation, but the values indicated from 
Figure 6 are in the range 10~ 4 to 10~ 5 , significantly lower than those used in Figure 5. The difference between spectra from 
models with these smaller values of r m and those from the Aizu (1973) formulation is less than 1% in the range of 0.05 eV 
to 100 keV, so that the observational effects of the leaky base will benegligible. Thus if we can indeed rely on small values of 
T m , current cold wall models such as those mentioned in Section 1 can continue to be employed for X-ray spectral fitting. 

The advantage of our treatment is that it avoids the infinities at the base of the flow inherent in all previous formulations 
in which the stationary-wall boundary condition base is assumed (e.g. Aizu 1973; Chevalier & Imamura; Wu 1994; Wu et al. 
1994). This is a step in the direction of more realistic formulations which allow the hydrodynamic accretion flow to match the 
hydrostatic white-dwarf atmosphere. As noted above the predicted spectral changes are small for small r m ; however a practical 
advantage is that it eliminates the numerical errors resulting from finite sampling of quantities which tend to infinity. In the 
previous treatments the fineness of the sampling of the base of the post-shock region in the numerical integrating scheme has 
an effect on the spectral slope: finer sampling will encounter a higher value of the density (which is tending to infinity) at the 
base, so that softer spectra will be generated. In our formulation, as the singularity at the base is removed, the effect of the 
sampling is insignificant even when small values of r m are assumed. This is consistent with the fact that emission from the 
boundary layer should be finite, as deduced from equation (45). 
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Figure 6. The values of T m (dotted lines) and (solid lines) as a function of the effective temperature of the white-dwarf atmosphere 
T c ff for 0.5- and 1.0-Mq white dwarfs (left and right panel respectively). The values of specfic accretion rates to calculate T h are 0.5, 1.0, 
5.0 and and 50.0 g cm s —1 (solid lines from left to right). Here, we have considered the mass-radius relation of white dwarfs given in 
Nauenberg (1972) and assumed that the temperature at the boundary T(0) is the same as T e g . The shock temperature T B is calculated 
by assuming a strong adiabatic-shock condition and a cosmic abundance. 
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A.l Convergence of the power series y") C n x n 

Consider a power series 
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Appendix A. 2 Accretion luminosity 

For the "leaky" accretion column with a specific mass flux a = — e~ T / Tm ), where /j, = (1 — e~ 1 / 4l " m )~ 1 , the accretion 
luminosity is 
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Hence, we have 
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